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Abstract 

An insulating optical lattice with double- well sites is considered. In the case of the 
unity filling factor, an effective Hamiltonian in the pseudospin representation is derived. 
A method is suggested for manipulating the properties of the system by varying the 
shape of the double-well potential. In particular, it is shown that the atomic imbalance 
can be varied at will and a kind of the Morse-alphabet sequences can be created. 



PACS: 03.75.Lm: 03.67.Hk: 03.67.Lx 



Keywords: Cold atoms, Double-well lattices, Order-disorder transitions, Quantum informa- 
tion processing 



1 



1 Introduction 



Physics of cold atomic gases, Bose [1-8] and Fermi [9,10], is a fastly developing field of 
research, both theoretically and experimentally. Cold atoms in optical lattices provide a 
versatile tool for various applications [11-14]. Recently, a novel type of optical lattices has 
been realized in experiments [15-17], the so-called double- well lattices. Each site of such a 
lattice is formed by a double-well potential. 

Generally, the theoretical description of the double-well lattices is essentially more com- 
plicated than that of the usual lattices. Therefore, in order to describe the properties of such 
lattices, it is useful to consider some particular types of them. 

In the present paper, we consider a particular kind of the double-well lattice, which 
is characterized by the following main features: atoms in the lattice are in an insulating 
state, and the atomic filling factor is equal to unity. Such a type of a lattice is of principal 
importance providing a convenient setup for realizing quantum information processing and 
quantum computing [11]. Another important property of the lattices, we shall be considering, 
is the existence of atomic interactions between different lattice sites. The principal goal of 
the present paper is to demonstrate that the atomic imbalance between the wells of each 
double well can be manipulated and that arbitrary sequences of the order-disorder transition 
can be generated. 



2 Insulating double-well lattice 

We start with the standard form of the energy Hamiltonian, expressed through the field 
operators in the Heisenberg representation. Since the system is assumed to be in the in- 
sulating state, the field operators can be expanded over localized orbitals, which yields the 
Hamiltonian 

H — ^EnC^jCnj + ^ X! ^jihhh 4c niji C li2j2 Cn 3h C n4j4 ; (1) 

nj {j} {n} 

in which c n j is a field operator labelled by the quantum index n and the site number j, 
related to a lattice vector a^ . The atoms can be either bosons or fermions, with the operator 
commutation relations 



where the commutator is assumed for bosons and anticommutator for fermions. The value 
E n represents the energy levels of an atom in a double-well located at a site of the lattice. 
The quantity ^Ij^hh^ * s a ma t rrx element of the interaction potential with respect to the 
localized orbitals labelled by the indices n and j. We assume that the interaction potential is 
sufficiently strong, so that the interactions of atoms between at least nearest-neighbor sites 
cannot be neglected. This can be easily achieved, for instance, with the long-range potentials, 
such that exist between polar molecules and between Rydberg atoms [18] or between the 
atoms with large magnetic moments [19]. A known example of the latter atoms is Cr that 
can be cooled to ultracold temperatures [20]. In the case of ions, long-range forces would be 
due to the Coulomb interaction. 
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The assumption that each lattice site contains just a single atom can be formalized by 
means of the unipolarity condition 

^ C nj C nj = 1 j c nj c nj = . (2) 

n 

This is the known condition, used earlier by Bogolubov [21] for treating ferromagnets. Keep- 
ing in mind low temperatures, we can consider only two lowest energy levels of a double-well 
potential, enumerated by n = 1,2, so that E\ < E 2 . As is known [22], the ground-state 
orbital, of an atom in a double-well potential, is symmetric with respect to spatial inversion, 
while that for the first excited state is antisymmetric. 

Retaining two lowest levels allows us to invoke the pseudospin representation by intro- 
ducing the pseudospin operators 



(c\jC2j - 4jClj) , S j =1 ( C l'C 2 j + <&jClj) ■ (3) 



qy ^ I 

3 ~ 2 \ Cl 3~*J -tj-'-J J i ~] 9 



To understand the meaning of these operators, one can introduce the left and right location 
operators, respectively, as 

1 1 

C jL = ^= (Clj + C2j) , C jR = -j= (dj - C 2j ) . (4) 

Then operators (3) become 

Sj = 2 ( c ]jL c jR + c )r c 3l) > S j = - 2 ( c ]l c jR ~ c )r c 3l) > 

S i = \ ( c UcjL - c) R c jR ) . (5) 

Note that this representation is valid for both statistics, so that atoms can be either bosons 
or fermions. The physical meaning of these operators is as follows: characterizes the 
tunneling intensity between the wells of a double-well potentials, Sj describes the Josephson 
current between the wells, and 5? defines the atomic imbalance between the wells. It is 
worth emphasizing that for a double well it is necessary to take into account at least two 
lowest energy levels, but not just one ground state, since only with two levels there exists 
tunneling between the wells. Actually, it is exactly because of the tunneling that the energy 
levels split into pairs [22]. 
Let us use the notation 

j-, _ Ei + E 2 
-^o = ^ , 



and 



where the sign plus or minus is for bosons or fermions, respectively. Also, let us denote the 
matrix elements 

A - Wt/1111 , t/2222 , OT/ 122l\ r> - Wt/1111 , t/2222 OT/ 122l\ 
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c* - \ iyr - V D > in - - 2 ^5 122 • ( 6 ) 

An important quantity is the tunneling frequency 

n = E 2 - E 1 + £ C iJ (7) 

characterizing the atomic tunneling between the wells of a double-well potential. 

Using the above conditions and notations, we reduce Hamiltonian (1) to the pseudospin 
form 

H = E N + \ ^ Aj - fi^SJ + £ lh 3 S?S* - ^hjS'S* . (8) 

The first two terms do not contain operators, so do not play role in what follows. The 
magnitude of can be comparable with Q, hence, it cannot be neglected. The value of 
the tunneling frequency Q can be varied in a wide range, depending on the shape of the 
double-well potential. To illustrate this, let us take the latter in the form 



' r \ 2 

W = v ( ^) 



To J 



+ V yz , 



where V yz is, say, a harmonic potential in the y- and z-directions. Then the tunneling 
frequency strongly depends on the parameter 



a = I = /^R ( 9 ) 

^2mr 2 V 2 V V ' 

in which a is the lattice spacing in the x-direction and Er = 7i 2 /2ma 2 is the recoil energy. 
By varying the interwell distance r or the potential depth Vb, the parameter a can be varied 
in a wide range between a < 1 and a ^> 1. Direct calculations [22] show that the tunneling 
frequency Q can also be made either small or large, such that 



n ~ 6\/ exp ^-2tt ^ ^LJ (a « 1) , 
/ a \ 2 / 3 /Ep\ 1/3 

^ 9,, »U (if) < a>>i >- < io > 

That is, f2 can be regulated by governing the shape of the double-well potential. And varying 
Q it is possible to regulate the system properties. 

The important system characteristics are the average tunneling intensity 



«■ = 4 £ < s* > , (n) 



N L 



average Josephson current 
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and the average well imbalance 



* = -£-E<s;>, (13) 

L 3 

where Nl is the number of lattice sites and < • > implies statistical averaging. These values 
depend on the system parameters, such as the dimensionless transverse interaction 

b = — — , (14) 
I+B ' V ; 

in which , ^ 

JVL i+3 JVL i+3 



and the dimensionless tunneling frequency 

us 7Tb- (15) 
3 Manipulating atomic imbalance 

The equations of motion for the average quantities (11) to (13) can be obtained by averaging 
the Heisenberg equations for the spin operators S". In this procedure, we assume that the 
system is at zero temperature, we employ the local- field approximation [23], taking into 
account particle interactions occurring in the local field of other particles. Measuring time 
in units of 1/(1 + B), we obtain 

dx dy 

— = (1 - b)yz - 7 2 (x - x t ) , ~fa = ( UJ ~ x ^ z ~ 7*(y ~~ Vt "> ' 
dz 

— = {bx - u)y - 7i(z - zt) , (16) 

where 71 and 72 are the attenuation parameters that are expressed through atomic interac- 
tions Iij and Bij similarly to the corresponding damping parameters in spin systems [24], 
and where the local fields are 

uj — bx 1 — 6 

x t = — ^ — , Vt = , z t = z , (17) 

with 



h= ^/{u - bx) 2 + (l-b) 



2 Z 2 



By their form, the local fields (17) correspond to the equilibrium solutions for averages (11) 
to (13). However, since here x = x(t) and z = z(t) depend on time, the local fields (17) 
describe the locally equilibrium state. 

Accomplishing the Lyapunov stability analysis for the system of equations (16), we find 
two fixed points. One is given by the equations 

x\ = u , y* = , z* = Vl-tu 2 (u < 1) , (18) 
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this fixed point being stable for u < 1, but unstable for to > 1. And the other fixed point is 

x* = l, y* 2 = 0, z* 2 =0 (w>l), (19) 

which is stable for uj > 1, but unstable for uj < 1. 

It is worth emphasizing that the fixed points (18) and (19) are stable, for the correspond- 
ing values of u>, only when ji and 72 are not zero. When the latter parameters are zero, 
the dynamical system (16) is structurally unstable. This stresses the necessity of taking 
into account the attenuation effects, without which there would be no correct description of 
dynamics. 

The stationary solution (18) characterizes the phase of the system with a nonzero well 
imbalance z\ > 0. Hence, this can be called the ordered phase. The stationary solution (19) 
describes a disordered phase, where the well imbalance is zero, z 2 = 0. In an equilibrium 
system, the transition between the ordered and disordered phases would correspond to a 
quantum phase transition, with the well imbalance playing the role of an order parameter. 
For a nonequilibrium system, the value uj = 1 is a bifurcation point, where a dynamical 
phase transition takes place. 

As has been discussed above, the magnitude of the tunneling frequency can be varied 
in a wide range by changing the shape of the double well. This means that we have a 
straightforward opportunity of regulating the state of the double-well lattice by reswitching 
the value of uj between that one corresponding to the ordered state and another related to 
the disordered state. That is, we can manipulate with the well imbalance, reswitching it, 
according to our will, between the zero and nonzero values. 

As an illustration of this remarkable feasibility of regulating the well imbalance, we solve 
numerically Eqs. (16), with a time-modulated tunneling frequency uj = u(t), which is varied 
according to the rule 

' u u < t < At x 
,.x ^2, At x <t < At x + At 2 

i uj!, At! + At 2 <t < Ati + At 2 + At 3 , 
^ 

where ui\ < 1 and uj 2 > 1- The time intervals At n can be chosen arbitrarily. For instance, 
we can take all of them being equal, as in Fig. 1, where At\ = At 2 = . . . = At n . Or we 
can organize a periodic sequence with two unequal time intervals Ati = At 3 = . . . = At 2n+ i 
and At 2 = At 4 = . . . = At 2n , as in Fig. 2. Finally, we can take arbitrary time intervals 
At n , as in Fig. 3, realizing the sequences, according to our will, in which information, as in 
the Morse alphabet, could be encoded. Such Morse-alphabet sequences for the pseudospin 
system, represent at ing a double-well lattice, are analogous to the punctuated sequences for 
superradiating spin systems [25]. 

In conclusion, we have derived a pseudospin representation for insulating double-well 
lattices, with the unity filling factor. Such a system can be either in an ordered or in a dis- 
ordered state, which are characterized by a nonzero or zero well imbalance, respectively. For 
correctly describing the system dynamics, one has to take account of relaxation effects, with- 
out which the related dynamical system is structurally unstable. By a temporal modulation 
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(20) 



of the lattice parameters, e.g., of the tunneling frequency, it is possible to regulate the well 
imbalance, organizing the Morse-alphabet sequences, which can be employed for quantum 
information processing. The possibility of regulating the well imbalance in a double-well 
lattice, by organizing arbitrary temporal sequences of ordered and disordered states, is the 
main result of the present paper. 
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Figure Captions 



Fig. 1. Population imbalance as a function of dimensionless time for the periodic 
reswitching of the tunneling between uj\ — 0.1 and uj 2 = 5, with equal time intervals Ati = 
At 2 = 10. The system parameters are 7i = 72 = 1 and b = 0.5. Initial conditions are 
x = 0.33, y = 0.8, z = 0.5. 

Fig. 2. Population imbalance as a function of dimensionless time for the periodic 
reswitching of the tunneling between lo\ = 0.1 and u} 2 = 11, with unequal time intervals 
Ati = 10 and At 2 = 25. Other parameters are as in Fig. I. 

Fig. 3. Population imbalance as a function of dimensionless time for a nonperiodic 
reswitching of the tunneling between uj\ = 0.1 and u 2 = 1.1, with varying time intervals At,-. 
Other parameters are as in Fig. 1. 
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Figure 1: Population imbalance as a function of dimensionless time for the periodic reswitch- 
ing of the tunneling between u\ — 0.1 and uj 2 = 5, with equal time intervals At! = At 2 = 10. 
The system parameters are 7i = 72 = 1 and b = 0.5. Initial conditions are x = 0.33, 
yo = 0.8, z = 0.5. 
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Figure 2: Population imbalance as a function of dimensionless time for the periodic reswitch- 
ing of the tunneling between cu 1 = 0.1 and u 2 = 1.1, with unequal time intervals Ati = 10 
and At 2 = 25. Other parameters are as in Fig. 1. 
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Figure 3: Population imbalance as a function of dimensionless time for a nonperiodic 
reswitching of the tunneling between u)\ = 0.1 and 002 = 1.1, with varying time intervals 
Atj. Other parameters are as in Fig. 1. 
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